Diffusive magnetotransport in a two— dimensional Rashba system 
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An analytical approach to calculation of the conductivity tensor, a, of a two dimensional (2D) 
electron system with Rashba spin-orbit interaction (SOI) in an orthogonal magnetic field is pro- 
posed. The electron momentum relaxation is assumed to be due to electron scattering by a random 
field of short-range impurities, which is taken into account in the Born approximation. An exact ex- 
pression for the one-particle Green function of an electron with Rasba SOI in an arbitrary magnetic 
field is suggested. This expression allows us to obtain analytical formulas for the density of states 
(DOS) and a in the self-consistent Born and ladder approximations, respectively, which hold true 
in a wide range of magnetic fields, from the weak (ui c t <C 1) up to the quantizing (lo c t > 1) ones. It 
is shown that in the ladder approximation the Rashba SOI has no effect at all on the conductivity 
magnitude in the whole range of classical (non quantizing) magnetic fields. The Shubnikov-de Haas 
(SdH) oscillation period is shown to be related to the total charge carrier concentration by the 
conventional formula, irrespective of the SOI magnitude. A simple equation defining the location 
of the SdH oscillation beating nodes is obtained. The results are in good agreement whith the 
experimental and recent numerical investigations. 

PACS numbers: 71.70.Ej, 72.15.Gd, 73.20.At, 73.21.-b 
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I. INTRODUCTION 

The growing interest in studying the spin-orbit in- 
teraction (SOI) in semiconductor two dimensional (2D) 
structures is mostly due to its potential application to the 
spin-based electronic devices^. There are two main types 
of SOI in the quantum well based on zinc-blende-lattice 
semiconductors. First, the Dresselhaus interaction^ that 
originates from the bulk inversion asymmetry (BIA); sec- 
ond, the Rashba interaction— induced by structural inver- 
sion asymmetry (SIA) of the confining field of a quantum 
well. Both of these interactions lead to the momentum- 
dependent spin splitting of the electron energy spectrum 
and to the formation of quantum states with the hard 
linked spatial and spin degrees of freedom of the elec- 
trons. They are responsible for many interesting effects in 
the transport phenomena like beatings in the Shubnikov- 
de Haas (SdH) oscillations^; weak antilocalizatio n 5 : 6 : 7 : 8 ; 
current-induced non-equilibrium spin polarizatiorAii; 
spin Hall eSeciM^, and so on. 

At present there are some sufficiently well developed 
theories of the kinetic and spin phenomena in 2D systems 
with SOI in zero or classical weak {uj c t <C 1) orthogonal 
magnetic fields. Here uj c — \e\B/mc is the cyclotron fre- 
quency, and t is the electron scattering time. As for the- 
oretical studies of the considered systems in strong, and 
especially in quantizing (u> c t > 1) magnetic fields, there 
is still no satisfactory analytical description of the kinetic 
phenomena even in the usual diffusive regime (without 
quantum corrections). The complex form of the eigen- 
spinors and energy spectrum of an electron in the pres- 
ence of SOI and a strong magnetic fields is the main cause 
of such a situation. Direct employment of this basis forces 
one to proceed almost right from the start to the numer- 
ical analysis of very cumbersome expressionsi^i^iSiSi. 



The strong magnetic field is however one of the most 
efficient tools for investigation of SOI^ and manipulation 
of the spin degrees of freedom in the semiconductor 2D- 
structures. Thus, a rather simple theoretical description 
of the kinetic phenomena in the 2D systems with SOI in a 
strong orthogonal magnetic field becomes a necessity. In 
the present work, we consider the problem of calculation 
of the longitudinal and Hall resistances of a 2D Rashba 
system in the ladder approximation assuming that the 
electron momentum relaxation is due to elastic scattering 
by short-range impurities which is taken into account in 
the Born approximation. 



We have found the exact relation between the one- 
particle Green function (GF) of the Rashba 2D-electron 
in an arbitrary orthogonal magnetic field and the well 
known GF of an "ideal" electron, that is an electron 
with the ideal value of the Zeeman coupling (g Q = 2) and 
without SOI. This allows one to obtain the analytical 
expressions for the density of states (DOS) in the self- 
consistent Born approximation (SCBA), and the conduc- 
tivity tensor a in the ladder approximation. The total 
DOS in SCBA is defined as the sum of the partial DOS's 
of two spin-split subbands. At the same time, conductiv- 
ity in ladder approximation looks as if the current were 
generated by charge carriers of one type with the total 
concentration n and mobility /z. These expressions hold 
good in a wide range, from the classically weak magnetic 
fields (lo c t <C 1) up to the quantizing ones (lo c t > 1). 
On the basis of these results, we perform the numerical 
analysis of the beatings of the SdH oscillations of the 
considered kinetic coefficients, as well as of their behav- 
ior in the classical magnetic fields region. The results are 
in good agreement with the experimental datai£*i*i, and 
with the results of a recent numerical investigation^ 7 ". 
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II. MODEL 

Let us consider a two-dimensional (||OAF) degener- 
ate gas of electrons with effective mass m, and effective 
Zeeman coupling g those move in an external orthogo- 
nal (B||OZ) magnetic field B = V x A in the presence 
of a random field U(r) due to pointlike impurities dis- 
tributed by Poisson law in the sample. We assume the 
Rashba interaction to be the dominant mechanism of the 
energy spin splitting in the absence of a magnetic field. 
This situation occurs, for example, in the narrow-gap 
semiconductor heterostructures, such as InAs/GaSb 4 , 
InGaAs/In Al AsiSiiS. The one-particle Hamiltonian of 
the considered system has the form 

7T 2 1 

5£ + Z7=- ha(er x 7r) n + -guj c a z +U(r) (1) 

2m 4 

(h = 1). Here tv = p eA/c = rov is the operator of 
the kinematic electron momentum; er = (a x , o~ y , a z ) is 
the vector formed by the Pauli spin matrices; a is the 
Rashba spin-orbit coupling; g is the effective Zeeman 
coupling (g-factor). 

In the gauge A = (0,Bx,0), the components of the 
eigenspinors of the Hamiltonian J{ (JIJ of a free (Z7(r) = 
0) Rashba electron are expressed through the Landau 
wave functions ip n x (r) depending on the Landau level 
number n = 0, 1, 2, . . . and the A-coordinate of the cy- 
clotron orbit center X — —k /muj~ 



f 



1 + C* 



C s>n ^ n -i,x(r) 



, a = (s, n, X) . 

(2a) 

The corresponding energy levels have the following form 



.[n + s^S 2 + 2 7 %] 



n = 0, s = +1 , 
n > 0, s = ±1 . 



_ (2b) 

Here C s . n = -/V2Ti/[sy/S 2 + 2 7 2 n - 8] is a normalizing 
coefficient; 5 = (g — 2)/4 is the relative deviation of the 
effective Zeeman coupling from its ideal value g — 2 (for 
definiteness, it is assumed that 5 < in these equations, 
but all the following results are valid for any sign of 5); 
and, finally, 7 = a^Jm/uj c is the dimensionless Rashba 
spin-orbit coupling. 

The quantum number s = ±1 describes the helicity 
of the Rashba electron eigenstate as in the absence of a 
magnetic fieldii. Indeed, it can be verified immediately 
that s = ±1 is the eigenvalue of the operator 



[acr x 7r + w c <5er] • n 

v /2ma 2 3i + uj 2 c 5 2 



(3) 



that is diagonal in the basis (|2a|l and approaches the 
helicity operator (er x p) • n/|p| as B -— > 0. Here n is the 
unit normal vector to the considered 2D system; 



'Kn 



7T 

2m 



1 



(4) 



is the Hamiltonian of the "ideal" (g = 2) electron in a 
magnetic field, which commutes with er ■ n, (er x it) ■ n, 
and with % QJ. 

In spite of this analogy with the B = case, we cannot 
say that the Rashba electron has in the states (|2a() the 
spin projection ±1/2 onto the direction a7r x n + w c 5n, 
because the components of the kinematic momentum op- 
erator 7r are not commuting motion integrals. Neverthe- 
less, this interpretation makes sense in the quasiclassical 
limit, when one can speak about the electron path in a 
magnetic field. Namely, the quantum number s = ±1 de- 
termines the value of the spin projection on the instant 
direction of an x n + oj c Sih that changes along the quasi- 
classical electron path. Thus, the spin configurations of 
the Rashba electron states form vortices in the AK-plane 
with center at the origin. 

The conductivity tensor a of the considered system 
has just one independent circularly polarized component 
a = a xx + io yx . In the one-electron approximation, it 
has the form 20 



G = -(- U 



8tt + 



E=E, 



(d E - d E 



4> 



AA 
EE' 



RR 
EE' 



dE 



E'=E 



(5) 



Here, $f|, = ( G x (E) V_ G Y (E')) is the current ver- 
tex operator; V± = V x ± iV y = v± ± 2iaa± are cir- 
cularly polarized components of the full velocity oper- 
ator [the corresponding components er are defined as 
<7 ± = (a x ± i<j y )/2], where the last term occurs due to 

SOI Q). G R( - A ~> (E) =l/(E-'K-U±iQ) is the resolvent 
(retarded (R) or advanced (A)) of the Hamiltonian Q, 
and angular brackets (. . .) denote the averaging over the 
random field U configurations. Finally, the symbol d E 
denotes the derivative with respect to energy E. 



III. ONE-ELECTRON GREEN FUNCTION 

By definition, the one-particle GF is the averaged re- 
solvent of the Hamiltonian © {G R{A \E)) = (1/(E - 
3i — U ± i0)). It is connected with the electron self- 
energy operator Y, R ( A > (E) by the relation (X = R, A) 



(G X {E)) = 



(G*(E)) (G^E)) 
(G*(E)) (G*{E)) 



1 



E-'K 



Z X (E) 
(6) 

The direct employment of the eigenspinors J2J) for cal- 
culation of @, or kinetic and thermodynamic proper- 
ties of the Rashba system in a strong magnetic field 
leads to very complicated expressions. One is forced 
almost from the first steps either to turn to numeri- 
cal calculations^*i£i^, or to make simplifying approxi- 
mations like the momentum-independent spin-splitting 
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energy 15 . This makes more difficult the interpretation 
of the results obtained in such a way, as well as the un- 
derstanding of the whole physical picture. But it turns 
out that the GF of the free (U = 0) Rashba system is 
expressed exactly through the GF of the " ideal" electron 
in a magnetic field. This opens up new possibilities for 
analytical studies of the considered system. Indeed, it is 
easy to check that the Hamiltonian of the free Rashba 
systems can be presented in the following form 



!K = 3-C + vy/2ma 2 Jf + to 2 c 5 2 



(7) 



Here v is the helicity operator defined in Eq. ©, JC is 
the Hamiltonian of the "ideal" electron 

The substitution of the Hamiltonian J7J) into the resol- 
vent G(E) = (E — gives, after some simple algebra, 
the following result (here and below, we drop superscripts 
R(A), if this does not lead to misunderstandings. Some- 
times, for brevity of notations, we shall not write explic- 
itly the energy arguments of the resolvents or GF's.) 

g = E - X + [q(7r x n) + u c 5n] ■ a 



(E + ma 2 — JC ) 



where 



n B = 2 v / 2ma 2 E + m 2 a 4 + u 2 S 2 = ^/fl 2 + Auj 2 c 5 2 . (9) 

The quantity Q B is equal to the magnetic field- 
dependent frequency of the spin precession of the elec- 
tron with energy E that is responsible for the Dyakonov 
— Perel spin relaxation mechanismSi; fl is the same fre- 
quency in the absence of a magnetic field. It should 
be noted that the same representation of the one 
electron GF can be also obtained for a system with 
the momentum-linear Dresselhaus SOI. For example, in 
the case of a [001]-grown quantum well based on the 
AjjjBy semiconductors, it is sufficient to replace tv — > 
7r = (iTy,ir x ) in the definition of the helicity operator J3J, 
change the sign before the Zeeman term (g = —2!) in 
the Hamiltonian of the "ideal" electron (0J, and, finally, 
to redefine the parameter 6 — ► 5 D = (g + 2)/4. 

The denominator of the right-hand side of Eq. JSJ) de- 
pends on the "ideal" electron Hamiltonian alone. Ex- 
panding this expression into the partial fractions, we ob- 
tain the desired representation of the one-electron GF of 
the free Rashba system 



G(E) 
1 



= E 

s=±l/2 



E 



n, 



--±1/2 



- 2s 



+ 4s [ma 2 — u c 8a z — a(iT x n) • cr] 

E + ma 2 

a(i x n) ■ a 



n, 



G(E + ma 2 + sfi B ) . 

(10) 



We use here the same notation (G) for the GF of the 
Rashba electron and for the GF of the "ideal" electron. 



However, this does not lead to confusion since the latter 
depends always on the energy arguments like E + ma 2 + 
sQ B etc. 

It is important that the same representation can be ob- 
tained for the averaged resolvent of the Rashba system in 
the SCBA. We restrict ourselves here to an approxima- 
tion in which the electron self-energy operator is diagonal 
in the spin space. Then, the SCBA equation for H X (E) 
has the following form 



t(E) = W(SpG{E)) = 



£ TT (£) 








(11) 



Here Sp denotes the trace only over the spatial degrees 
of freedom; W = rijU 2 , where n 1 is the impurity concen- 
tration, £/q is the magnitude of the pointlike potential of 
an isolated impurity. Therefore, it is sufficient to make 
everywhere in Eq. IjlOjl the following substitutions 



E -» E - E e (£) 



9^c 



gu c + 4Z (E) (12) 



to obtain the desired representations for the averaged 
GF's in the SCBA. Here S e(o) (E) = [E n (E) ± (E)] /2 
are the even and odd parts of the electron self-energy. 
The first (S e = A e ± i/2r e ) describes the perturbation 
(shift A e and broadening l/r e ) of the one-electron energy 
levels by a random field. The real part of S G = A ±t/2r 
defines the renormalization of the Zeeman coupling i|12fl , 
while its imaginary part oc 1 jr makes a contribution to 
the overall broadening of the one-electron energy levels. 
As a result, we obtain a expression like Eq. IjlOj l for the 
averaged GF, where 



1 



G R{ - A \E + ma 2 + sn B ) = 

E + ma 2 + sfl B - J£ ± — 

(13) 

is the averaged retarded (advanced) GF of the "ideal" 
electron, and 



1 1 • fr^R r^A\ 4ma 2 \ 1 Au) c 5 1 

t„ r e ' \ n B j r e n B t 

(14) 

are the disorder-modified frequency of the spin preces- 
sion 10 and the inverse life time of an electron in the s-th 
spin-split subband. As usual, we do not take explicitly 
into consideration in Ijl3jl the one-electron energy levels 
shift A e that is absorbed by the normalization condition, 
but we mean here that the odd shift A Q is included in 
the definition of the effective ^-factor in accordance with 
(j!2jl . The explicit allowance for the Zeeman coupling 
renormalization is particularly important in the SdH os- 
cillations regime. 
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IV. DENSITY OF STATES AND SELF— ENERGY 

We first consider the calculation of the DOS n(E) = 
Im(Tr G A (E)) /it using the above-obtained expression for 
the one-particle GF (|10|l . Here, the symbol Tr denotes 
the trace over the spatial and spin degrees of freedom. 
For the sake of simplicity, we shall deal with the case of 
large filling numbers (E 3> w c ). Calculating the trace of 
resolvent (ITUf) over the spatial and spin degrees of free- 
dom, wc obtain the following expression for the DOS 



n(E) = V ^n^\E + ma 2 + s(Q b ±lu c 

s=±l/2 

= E ^ 0) (£). 

t-^ m 

s=±l/2 



(15) 



Here, we take into account that the DOS of a spinless 
electron in an orthogonal magnetic field n^ ' (E) satisfies 
n^(E) = n^(E ± lu c ) at large filling factors (E > uj c ). 
The sign before uj c is chosen in such a way as to ensure 
the right-hand limit s(VL B ±us c ) — + ±sgu} c /2, as the spin- 
orbit coupling approaches zero. The effective mass m s in 
the s-th subband is defined as 



m 1 



Araa" 



md E (E + sVt B ). (16) 



In the considered case this expression coincides with the 
usual definition of the transport and cyclotron effective 
masses in the isotropic nonparabolic band^. 

In full accordance with the two-subband model, the 
DOS in Eq. I|15|) is presented as a sum of partial con- 
tributions. Using this expression for the DOS, we can 
obtain the analytical form of the equation for the elec- 
tron concentration n = J F n(E)dE that is the normal- 
ization condition for the Fermi level determination. So, 
a more correct expression for the DOS is needed, which 
adequately describes its behavior not only in the vicin- 
ity of the Fermi level E F (> 0), but also near the lower 
boundary of the spectrum. For example, at B = we 
have 



ma 



n(E) 



yj2mE + m 2 a 2 
1, 



1 9 

--ma z < E < , 
E > 0. 



(17) 

In the energy interval — ma 2 /2 < E < 0, the DOS is 
formed by the states of the lower spin-split subband and 
has the typical one-dimensional behavior. Integrating 
(f 1 71) between —ma 2 /2 and E F , we obtain 



(E p + ma 2 ) = — E , E F > . 



(18) 



Thus, the energy E Q = E F 



corresponds to the 



Fermi level in the absence of SOI. Notice that the partial 
electron concentrations n s = m(E + sfl B ) /2tt depend 
nonlinearly on the Fermi energy, in contrast to n (JTSJ. 
The correction to the Fermi energy in (|18fl comes from the 



low-energy tail of the DOS (|17|l . Of course, the difference 
between E and E F is small for weak SOI (ma 2 <C E F ). 
However, as shown below, it is of crucial importance to 
take it into account for a correct interpretation of the 
spin-orbit interaction effect on both the conductivity in 
the absence of a magnetic field, and the SdH oscillations. 

The representation i|15|) allows one to obtain a simple 
analytical expression for the DOS that holds good up to 
the quantizing fields region (uj c t > 1). Indeed, the DOS 
of a spinless electron in the large filling factors region 
(E ^> uj c ) has the form 



sinh ■ 



in 



cosh h cos 2-7T — 



(19) 



Inserting Eq. (|19|) into Eq. I|15|) . we obtain for the oscil- 
lating part of the DOS the following expression 



» . . 2m 
An(E F ) = exp 

7T 



7T 



En n 

COS27T— - C0S7T — — 



2ma 



e, 



■ sin 2n — - sin7r 



(20) 



that is valid in the magnetic fields region under consid- 
eration. The second term in Eq. 1)2(1 ll appears due to the 
difference between the effective masses m s <|16l) . 

It follows from 1)20(1 that energy E defines the main pe- 
riod of the SdH oscillations. In other words, in the large 
filling factors region (E F ^> lu c ) the SdH period is deter- 
mined by the total electron concentration (see Eq. H18|) ~). 
redardless of the spin-orbit interaction magnitude. On 
the other hand, the period of the SdH oscillation beatings 
(|20|1 (see Fig. 1(a)) defines the spin precession frequency 
@ which depends on the magnetic field even in the ab- 
sence of Zeeman splitting (g = 0). 

In the case of weak SOI (ft -^i E), the oscillations of 
the DOS are determined completely by the first term in 
Eq. (|2U|I . Then, the location of the fc-th node of beatings 
is determined by the condition 



B k = 



2mc 



\e\ v /(2fc + l) 2 -(< 7 -2)2 



(21) 



This limit was considered in Ref. |2ji Unlike the re- 
sults of that work, the above-obtained equations still 
stand in the case of strong SOI, where it is important 
to take account of the difference E Q and E F for correctly 
defining the SdG oscillation period. In addition, we have 
taken into account the Zeeman splitting of the electron 
spectrum that allows to describe more correctly the os- 
cillation pattern. For example, the Eq. 1)21(1 allows us to 
determine both the spin-orbit, a, and Zeeman, g, cou- 
plings by the measured locations of two different nodes 
(see the upper curve in Fig. 1(a)). On the other hand, 
the spin precession frequency Q B approaches |<5|o; c as the 
magnetic field B increases. Therefore, in this case a grad- 
ual transition from the beatings of the SdH oscillations 



5 




partial DOS's with opposite spin projections onto the 
OZ-axis 



COT 



FIG. 1: Plots of the SdH oscillations of the total DOS (up- 
per panel) and of the difference of the partial DOS's (bot- 
tom panel) of the 2D Rashba system at fixed g — 2.8, 
and k F l = 35.0, and different fir = 3.0; 1.5; 0.75 (up to 
down). The arrows point the nodes location with their num- 
bers k = 2, 3, 4, 5 that are calculated with Eq. 1211 . 




FIG. 2: Plots of the SdH oscillations of the inverse life time 
of one-electron states in the s-th spin-splitted subband at 
different values of Zeeman factor g — 1.8; 1.0; 0.2 (up to 
down), and fixed values of k F l = 35.0, and fir = 1.5. 



8n{E) = n n (E)-n ll (E) 



4uJ 



sn^(E) (22) 



= ±1/2 



This quantity is proportional to the derivative of the 
transverse spin magnetization with respect to energy E 
and, therefore, it enters in the definition of the effective 
concentrations of current carriers in the dissipative part 
of the 2D Rashba system conductivity in an orthogonal 
magnetic field (see the next section). 

Evidently, the difference of the partial DOS's (12211 van- 
ishes in the region of classical magnetic fields (oj c t C 1), 
but it plays an important role in the SdH oscillations 
regime. In the case of large filling factors, the oscillating 
behavior of this quantity is described by the following 
expression 



5n(E f 



2m 2lo„5 



■ exp 



• „ E o ■ ■ 
sm27r — siH7T- 



(23) 

Unlike the total DOS l|2U [) . this expression contains just 
one oscillating term, because Sn(E) does not depend on 
the effective masses m s (|lfci|> . Indeed, the difference of 
the partial DOS's 8n(E) is non zero, which is entirely 
due to the spin degrees of freedom of the electrons. The 
typical SdH oscillation patterns of Sn(E) are depicted in 
Fig. 1(b). 

Now, let us turn to the discussion of the electron life 
time t s in the s-th spin-split subband which is defined, 
according to Eq. Ijl4(l , by the imaginary parts of the even 
and odd self-energies S , v In other words, the total life 
time of the one-electron states r s is determined by the 
sum of the weighted relaxation rates of the orbital and 
spin degrees of freedom. The first term in this expression 
is proportional to the above-considered total DOS, hence 
its magnetic-field dependence coincides up to the scale 
factor with the patterns shown in Fig. 1(a). Of particular 
interest is the last term in Eq. (|14fl stemming from the 
Zeeman coupling renormalization. It is proportional to 
the difference of the partial DOS's (1221 and, therefore, 
plays an important role in the SdH oscillation regime, as 
shown in Fig. 2. Notice that the beatings of the SdH 
oscillations are supressed with the increase of the rela- 
tive magnitude of the second term in Eq. (|14|) . Indeed, 
Eq. H21(l determines the location of the beating loops of 
the oscillation instead of the nodes. Thus, the broaden- 
ing of the Zeeman levels leads to observable supression 
of the beatings of the electron life time t s oscillations. 



to the familiar Zeeman splitting of the oscillating peaks 
should be observed. The beginning of this transition can 
be seen on the lower curve in Fig. 1(a). 

Another important characteristic of the one-electron 
states of the 2D Rashba system is the difference of the 



V. CONDUCTIVITY 

The general expression for the conductivity (JSJ con- 
sists of two different terms. The first of them describes 
the contribution of the electrons at the Fermi level, the 
second one contains the contributions of all filled states 
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below the Fermi level. We begin the calculation of the 
conductivity with the last term of a 11 . First of all, it 
is pure imaginary and, therefore, makes a contribution in 
the Hall conductivity alone. Stfeda™ was first to show 
that, for spinless electrons, this part of the conductivity 
is equal to 



i\e\c 



dn 
~dB 



(24) 



interpretation can be found in survejiSi. 

This result is extended immediately to the electron sys- 
tems with SOI. Following Stfeda's argument, it can be 
shown that the part a 11 of the 2D Rashba system con- 
ductivity is expressed as 



a 11 =i\e\c 



dn\ (dM v 



(25) 



where n is the electron concentration. It should be 
pointed out that Eq. (|24|) is exact, and with the thermo- 
dynamic Maxwell relation a 11 can be expressed through 
(dM/dE) B , where M is the orbital magnetization of the 
electron gas. Detailed discussion of a 11 and its physical 



where M p is the spin magnetization of the electron gas. It 
follows that in the general case a 11 is determined only by 
the diamagnetic part of the electron gas magnetization. 
By direct differentiation of the electron concentration n 
with respect to the magnetic field induction B we obtain 



J 



dn 
dB 



n 
B 



Tl ^ A _ dE = n_ 
dB v 1 2ni B 



Tr 



d'K dG 1 



dZ A dG / 



dY, A dG / 



dB dE dB dE dE dB 



h.c 



where the symbol h.c. denotes the Hermitian conjugate 
the terms. Considering that in the SCBA S = W Tr G/2, 
it is easily seen that the expression in the square brackets 
is the total derivative with respect to energy. As a result, 
integration by E in 126|) is performed explicitly, and after 
some simple algebra the expression for a 11 takes on the 
form 



where 



(27) 



(28) 



Eq. H27|l is a generalization of the known SCBA expres- 
sion for a 11 in the case of spinless electrons^. Particu- 
larly, the quantity n ± (|28|l is the counterpart of the famil- 
iar parameter n ± = EnS°\E) that stands for the current 
carrier concentration in the dissipative part of the con- 
ductivity tensor of spinless 2D-electrons in a magnetic 
field in the SCBA 20 . In the classical (i.e. nonquantizing, 
lu c t <C 1) magnetic fields region, it is equal to the total 
electron concentration n. 

Now, we turn to the first term in the conductivity JSJ. 
It is quite easy to show, by identical transformations, 
that 



8tt 



TrV + [*jA+*™] 



Attiti 



Tr (G A + G 1 



(29) 



The main contribution to the dissipative part of the con- 
ductivity is proportional to the current vertex 'I'^ in 
Eq. (JSJ. If we accept the SCBA {TTJ for the electron 




V 



(a) 

V_ 



V 



(b) 



V 



FIG. 3: (a) The diagrams depicting the conductivity in the 
ladder approximation; (b) The impurity ladder series in the 
particle-hole channel (diffuson). 



self-energy £, we must evaluate this part of conductiv- 
ity in the ladder approximation in order to satisfy the 
particle conservation law (see. Fig. 3). 

Replacing, as a first approximation, TrVj_$^ — > 
we obtain the "bare" conductivity 



TrV + G R V„G A 



= — Tr 

47T 



V + G R V_G A + — (G / 
m 



G 



, (30) 



depicted by the first diagram in Fig. 3(a). In what fol- 
lows, for simplicity, we neglect everywhere the odd part 
of the electron self-energy E G l|14f> . In this approximation 
one can obtain the following equations 



7 



-' + — ] G v_G = r r;- 1 g't. + ->h, 

^)g a V+G r = G a V„ - WO 1 ' - 2/n 



r 



<r_G - G J V 



G R a_G J 



G a, - a,G G a,G 

T 



(31) 



Using these relations, we perform a series of transfor- 
mations of the "bare" conductivity l|3Tj|) . neglecting the 
terms being small in parameter 1/E f t and lying beyond 
the accuracy of the ladder approximation. Omitting in- 
termediate transformations of purely technical character, 
the final result can be written as 



a lare 



e 2 r„ 



1 



m 1 



2ma 2 n F 
1 — iu} n T< 



■{1-P) 



(32) 



where n F = Tr(G A - G R )/4m is the DOS at the Fermi 
level per spin; 



P 



W Tr a + G a_G J 



(33) 



The ladder correction to the conductivity shown by the 
second diagram in Fig. 3(a) is presented by the following 
analytical expression 



Aa/ ad = — tr 



4-7T 



Sp (G A V + G R ) ■ f> • Sp (G R V_G J 



(34) 

The "bare" current vertices Sp {G A ^V ± G Ri - A ^) in- 
volved in H34|l can be expressed by the parameter P Ij33(l 
and the relaxation time r e , using relations l]3ip. In the 
main order in the small parameter 1 / E f t they have the 
form 

Sp {G A ^V ± G R ^) = Tt 47ran f T e (i _ p) ^ 

1 llO c T e 

The diffuson D depicted in Fig. 3(b) by a ladder se- 
ries in the particle-hole channel is a 4 x 4-matrix in the 
representation of the total spin of the electron-hole pair. 
The Bcte-Salpeter equation for T) (see the second row in 
Fig. 3(b)) has the same structure. However, as seen from 
(|34|l and (|35|l . the contribution to Acr^ comes from the 

scalar quantity CD = trcr, Da_ being the projection of 
the diffuson on the triplet state 1 1 , — 1) of the electron- 
hole pair. Projecting the Bete-Salpeter equation on this 
state, we obtain a closed scalar equation for D with the 
solution of the form 



D 



W 
1 - P ' 



(36) 



where P is defined in 133|) . 

Substitution of (|3*6*j) and (|S5|l into (|^4*jl yields the ex- 
pression for the ladder correction to the conductivity 
which exactly cancels the second term in the square 
brackets in i|32|l . As a result, the final expression for 



the conductivity of the Rashba system in a transverse 
magnetic field in the ladder approximation takes on the 
form 



B 



1 - iu c T e 



(37) 



Eq. ipT7|) looks as if the current were generated by 
charge carriers of one type with mobility /i — \e\T e /m 
and concentration n. This would be expected, because 
the conductivity tensor in the absence of a magnetic field 
is diagonal in the original spin space (a^ = a ^ = 0) by 
virtue of the momentum parity of the GF's, and the full 
conductivity is equal to a = + c^™- In the classical 
(nonquantizing) magnetic fields this property is retained, 
since the difference of the partial DOS's with opposite 
spin projections onto the OZ-axis (|22|l is equal to zero 
in this region (see Fig. 1). 

Of course, the simple structure of Eq. (|37|1 for the con- 
ductivity breaks down in the region of sufficiently strong 
magnetic fields, uj c t 1, where the contribution l/r 
to the inverse life time of the one-electron states <|14[) 
cannot already be neglected. As shown in the previous 
section, this may lead to in flattening of the SdH oscil- 
lation beatings (see Fig. 2). Therefore we believe the 
approximation used in this section 1 jr a = to be correct 
in the magnetic field region where well-defined beating 
nodes of magnetooscillation are observed. 



VI. RESULTS AND DISCUSSION 

First of all, let us summarize briefly the main results 
obtained in this work. We have shown that the eigen- 
states of the 2D Rashba electron in an orthogonal mag- 
netic field are characterized by a special motion integral 
© that generalizes the notion of helicity^. Using this 
fact, we have found the relation IjlUI) between the GF's of 
the 2D Rashba electron and the "ideal" one that holds 
good for arbitrary orthogonal magnetic fields as well as 
for the strong spin-orbit coupling. With the help of this 
relation, we have obtained, in contrast to Refs. ll4Tll5lll^ . 
the analytical expressions for the DOS in SCBA (|15() and 
for the magnetoconductivity in the ladder approximation 
(|37|l of the 2D Rashba system that are valid in a wide 
range from the classical magnetic fields up to the quantiz- 
ing ones (ijJ c t ^1)- The spin-orbit as well as the Zccman 
splitting of the electron energy are properly a llowed for 
in these expressions, unlike the results of Refs. I23l2 ll In 
particular, we have obtained a simple expression (|21|) for 
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the nodes locations of the SdH oscillations beatings. We 
have shown that the competition of the relaxation rates 
of the orbital and spin degrees of freedom in the total in- 
verse life time 1 /r s of the one-electron states in the s-th 
subband leads to the partial supression of beatings of the 
1/t s SdH oscillations. 

We start the discussion of the results with the conduc- 
tivity in the classical magnetic fields region (uj c t c < 1). 
In this case, it follows immediately from Eqs. I|28[l . and 
(|37|l that the conductivity of a 2D Rashba system takes 
the usual Drudc Boltzmann form 



D 



where 



1 - iu c T e 



-(E F + ma 2 ) 



(38) 



(39) 



is the Drudian conductivity in absence of the magnetic 
field; r e = 1/ mW is the life time of an one-electron state 
at B = 0. It immediately follows that in the ladder 
approximation the classical magnetoresistance of a 2D 
Rashba system is zero, p(B) = p D = l/o~ D , and the Hall 
coefficient R H = — l/\e\nc. 

Thus, in the ladder approximation the Rashba spin- 
orbit interaction has no effect at all on the conductivity 
magnitude everywhere over the region of classical mag- 
netic fields, uJ c T e <C 1 (including the case of B = 0). 
Note that the first relation for the Drudian conductiv- 
ity (|39|l formally coincides with that obtained in Ref. I2H 
However, as mentioned above, the correction ma 2 to the 
Fermi energy does not lead to an observable change in a, 
since it is absorbed by the normalization condition (|18H . 

Let us proceed now to the discussion of the magneto- 
transport in the 2D Rashba system in the the large filling 
factors (E ^> uj c ) region, where the SCBA and the lad- 
der approximation are applicable to the description of the 
one-electron states and kinetic phenomena, respectively. 
As usual, we extract in the linear approximation the oscil- 
lating parts of the conductivity that enter through DOS 
into the effective concentration <|28[l and mobility p. 
As in the ladder approximation the conductivity (|37|l has 
the form characteristic of a conductor with one type of 
charge carriers, we can immediately use the expression 
obtained in Ref. [25j for the oscillating parts of the longi- 
tudinal resistance p and the Hall coefficient R H 



Ap(B) _An(E F ) 



Po 
R° 



1 An{E F ) 
< 2 t 2 „(f>) 



(40a) 
(40b) 



Here, p = 1/<J D and R° H = —l/\e\nc are the resistance 
(see Eq. (j2HJ)) and Hall coefficient in zero magnetic field, 
respectively. Thus, in the linear approximation the SdH 
oscillations p and R H are entirely determined by the 
magnetic-field dependence of the DOS, An(E F ), the first 



harmonics of which are of the form (|20|l . From this we 
can draw two conclusions that are of great importance for 
an adequate interpretation of the SdH oscillation pattern 
in systems with spin-orbit interaction. 

First, the period of the SdH oscillations of p and R H is 
defined by energy E = E F + ma 2 , and not by the Fermi 
energy E F , as was stated in Refs. I2I2I From this and 
the normalization condition l|18fl it follows that the SdH 
oscillation period is related to the total charge carriers 
concentration n by the well known formula 



B 



(41) 



that holds true irrespective of the magnitude of the spin- 
orbit interaction constant. 

Second, the SdH oscillation beating period is defined 
by the magnetic-field-dependent spin precession fre- 
quency fl B = y/fl 2 + (g - 2) 2 w 2 (see. Eq. ©). From 
this it follows a simple equation for the beating nodes 
location B k (|21(l that allows the constants a and g to be 
found from the measured B k values. This is illustrated 
Fig. 4 which presents the results of least-squares fitting of 
Eq. lf?T|l in variables B~ 2 , (2fc + 1) 2 to the measured^ lo- 
cations of the SdH oscillation beating nodes of the longi- 
tudinal resistance of a In 2 .Ga 1 _ x As/In 52 A1 48 As-type 
heterostructure. The slope of the fitted straight line and 
the point of its intersection with the ordinate axis yield 
the values f2 = 2.47 meV and g — 4.25 for the spin pre- 
cession frequency in the absence of magnetic field and the 
^-factor, respectively. This values are in good agreement 
with the results £1 = 2.46 meV and g — 4.4±0.2 obtained 
in Ref. [H 

It should be stressed that fi B = O w 2k F a only at 
g = 2. Only in this case the beating nodes locations 
are strictly periodic in the reverse magnetic field and are 
described by the condition 2k F a = (k + l/2)w c obtained 
in Refs. I2II2I 

Using Eq. I|37l) . we have performed numeric analysis 
of the SdH oscillations of the longitudinal resistance p 
(Fig. 5) and the Hall coefficient R H (Fig. 6) for param- 
eters (g = 4, E = 108.93 meV, O = 5.13 meV, and 
g = 3.8, E = 98.85 meV, fl = 5.59 meV) correspond- 



ing to the gate voltages V — 0V and V 



-0.3V 



for the In Ga As /In Al As heterostructure investigated in 
Ref. The parameters g and O were calculated from 
the positions of two successive nodes B\, and B2, using 
Eq. (|21|l . while Eq was adjusted to the SdH oscillation 
period. At V g — V, the calculated value f2 = 5.13 meV 
is close to the result fl = 5.4 meV of Ref. 18, whereas 
at V g = —0.3 V, the calculated value = 5.59 meV is in 
excellent agreement with the value = 5.6 meV obtaned 
in Ref. The results of p calculation are compared in 
Fig. 5 to the experimental curves from Ref. IT3. It can 
be seen that the theoretical results reproduce well the 
period and beating nodes location of the measured mag- 
netoresistance oscillations. Some difference in oscillation 
amplitude is due to the fact that the temperature smear- 
ing of the Fermi level was not taking into account in our 
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50 




-10 1 1 1 1 1 1 1 ' ' 1 

20 40 60 80 100 120 140 160 180 

(2k+1) 2 

FIG. 4: Measured locations, B k , of the SdH oscillation beat- 
ing nodes as compared to the behavior predicted by Eq. 12H . 
The points correspond to the values Bi = 0.873 T, B 2 — 
0.460 T, B 3 = 0.291 T, B 4 = 0.227 T, B 5 = 0.183 T, B 6 = 
0.153 T, measured for the Ino.65Gao.35 As/In .52 Al .4gAs- 
heterostructura— . The straight line is the result of least- 
squares fitting of Eq. l lSfl . 




2 [ 1 1 , 1 1 1 , 1 1 , , 1 1 . 1 . . 1 

0.5 1 1.5 2 2.5 

B (Tesla) 

FIG. 5: The theoretical curves of the p magnetooscil- 
lation as compared with the measured— results for the 
InGa As/In AlAs-heterostructure at gate voltages V g — V 
and V g — —0.3 V. The experimental data are denoted by the 
thick solid line and the theoretical ones by thin solid line. 



analysis for simplicity. The negative magnetoresistance 
observed in Ref. [lilies outside the ladder approximation. 

In conclusion it should be stressed once more that the 
results of this work have been obtained with rigorous 
account taken of both the spin-orbit and Zeeman split- 
ting of the energy levels. Up till now this was made 
by using numeric analysis only 14-16,17 . Our results re- 
produce quantitatively all the peculiarities of the mag- 
netooscillation curves obtained in such a way for models 
with Rashba spin-orbit interaction 1 ^ 7 -, except for the 
anomalously large positive magnetoresistance obtained 
111 Ref. B Our conclusion about the absence of pos- 
itive magnetoresistance is in drastic contradiction with 



0.5 








0.5 1 1.5 2 2.5 

B[T] 



FIG. 6: Plots of the SdH oscillations of the Hall coefficient 
of the 2D Rashba system calculated for the parameters corre- 
sponding to the InGa As/In Al As-heterostructure— for gate 
voltages V g — V and V g = —0.3 V (from the bottom up- 
wards) . 



Ref. Ilil but agrees with the results of analyticalS&Sl and 
recent numerical^ studies. 
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